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ABSTRACT 

This paper presents the first optical variabihty study of the Westerlund 1 
super star cluster in search of massive eclipsing binary systems. A total of 129 
new variable stars have been identified, including the discovery of 4 eclipsing 
binaries that are cluster members, 1 additional candidate, 8 field binaries, 19 
field 6 Scuti stars, 3 field W UMa eclipsing binaries, 13 other periodic variables 
and 81 long period or non-periodic variables. These include the known luminous 
blue variable, the B[e] star, 11 Wolf-Rayet stars, several supergiants, and other 
reddened stars that are likely members of Westerlund 1. The bright X-ray source 
corresponding to the Wolf-Rayet star WR77o (B) is found to be a 3.51 day 
eclipsing binary. The discovery of a reddened detached eclipsing binary system 
implies the first identification of main-sequence stars in Westerlund 1. 

Subject headings: binaries: eclipsing, stars: Wolf-Rayet, stars: variables: other, 
open clusters and associations: individual (Westerlund 1) 



1. Introduction 



Westerlund 1 is a unique laboratory for studying the elusive evolution of massive stars. 
It has recently emerged as the nearest known super star cluster, i.e., a glob ular cluster 
precursor, typically found in starburst galaxies. Its discovery dates back to IWesterlund 
(Il96ll ). however, the high amou nt of reddening ( Ay ~ 12 mag) hi ndered spectroscopic 



observ ation s for many 



y ears (see 



Westerlundl 119871 ). Spectroscopy by IClark fc Negueruela 



(I2OO2I ) and IClark et al.l (120051 ) revealed a massive stellar population of post-main sequence 
objects and determined it to be the n iost m assive compact young cluster known in the Local 
Group. iMengel fc Tacconi-GarmanI (120071) measured a dynamical mass of logM = 4.8, in 
agreement with the mass inferred by iGlark et al.l (120051 ) . Westerlund 1 therefore weighs more 
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than the R136 or Arches clusters, however, with an inferred age of 4.5-5 Myr fICrowther et al. 
20061 ). it is more similar in age to the Quintuplet and Center cluster. Furthermore, it is less 
obscured by dust than the Galactic Center clusters, which allows for observations in the 
optical and is ten time s nearer tha n R136 , which decreases crowding and blending effects 
by the same amount. IClark et al.l (120051 ) obtained spectroscopy of 53 stars, which were 
all classified as post main-sequence objects, including blue supergiants, red supergiants, 
yellow hypergiants, a luminous blue variable and 14 ^ yolf-Rayet (WR) star s. Currently, 24 



WR stars have been identified in Westerlund 1 (see iCrowther et al.l |2006| . and references 
there in), concentrating 8% of the known galactic WR stars in one cluster (Ivan der Hucht 
20061). Chandra o bservations have revealed the presence of a magnetar in Westerlund 1 
(IMuno et al.l l2006l ) . which for the first time places a firm lower limit on the mass of a 
neutron star progenitor. In addition, 12 WR stars and a B[e] star were found to be X-ray 
sources (ISkinner et al.ll2006l ). 



The motivation behind this paper is twofold: to discover eclipsing binaries containing 
massive stars (> 5OM0) and to provide time domain information for this remarkable cluster. 
Eclipsing binaries are extremely powerful tools that can be used to measure the most mas- 
sive stars, probe the upper stellar mass limit and provide accurate masses and radii for the 
most massive stars in a variety of environments and at a range of metallicities. Double-lined 
spectroscopic binary systems exhibiting eclipses in their light curves provide accurate geo- 
metric measurements of the fundamental parameters of their component stars. Specifically, 
the light curve provides the orbital period, inclination, eccentricity, the fractional radii and 
flux ratio of the two stars. The radial velocity semi-amplitudes determine the mass ratio; the 
individual masses can be solved using Kepler's third law. Furthermore, by fitting synthetic 
spectra to the observed ones, one can infer the effective temperatur es of the stars, solve for 
their luminosities and derive the distance (e.g. iBonanos et al.l 120061 ). 



Until recently, the most massive stars measured in eclipsi ng binaries were R136-38 in 
the Large Magellanic Cloud (56.9 ± 0.6 M^. iMassev et al.l [20oi) and WR 22 (55.3 ±7.3 M©, 
Rauw et al.l Il996l : ISchweickhardt et al.l Il999l ) , an evolved star in our Galaxy. The current 



heavy weight champion is a galactic WR binary, WR 20a (iRauw et al.ll2004l : iBonanos et al. 



20041 ). in the young compact cluster Westerlund 2, with component masses 83.0 ± 5.0 Mq 
and 82.0 ± 5.0 Mq, making this the most massive binary known with an accurate mass 
determination. Such systems are of particul ar interest, since massive binaries might be 
progenitors of gamma-ray bursts ( GRBs, e.g. Frver et al.| 2007 ). especially in the case of 
Population III, metal- free stars (see iBromm fc Loebll2006l ). They also provide constraints 
for untested massive star formation and evolution models, stellar atmosphere and wind 
models. Acc urate parameters of massive stars provide insights into the frequency of "twin" 
binaries (see iPinsonneault fc Stanekll2006l : iKrumholz &: ThompsonI 120061 ) . the progenitors of 
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X-ray binaries, core-collapse supernovae, the connection between supernovae and GRBs, and 
Population III stars, by studying any observed trends with metallicity. 

A systematic search for the most massive stars (> 50Mq) in eclipsing binaries is currently 
underway. Analogs of WR 20a, if not more massive binaries, are bound to exis t in the young 



mass ive clusters at the center of the Galaxy (Center, Arches, Quintuplet; e.g. iPeeples et al. 



20071 ). in nearby super star clusters (e.g. Westerlund 1, NGC 3603 and R136), in Local Group 



galaxies (e.g. LMC, SMC, M31, M3 3) and beyond (e.g. M81 , M83, NGC 2403). The large 



binary fraction (> 0.6) measured by iKobulnicky et al.l (l2006l ) among early-type stars in the 



Cygnus 0B2 association implies that a significant number of massive short-period eclipsing 
binaries should exist in massive clusters. The first step in measuring accurate masses of 
the most massive stars involves discovering eclipsing binaries in massive clusters and nearby 
galaxies. The brightest of these are expected to contain very massive stars, and follow-up 
spectroscopy provides fundamental parameters for the component stars. 

This paper presents the first variability study of the massive Westerlund 1 cluster, in 
search of massive stars in eclipsing binaries. It is organized as follows: §2 describes the 
observations, §3 the data reduction, calibration and astrometry, §4 the variable star catalog, 
§5 the color magnitude diagram and §6 the summary. 



2. Observations 



Photometry of Westerlund 1 was obtained with the Direct CCD on the 1 m Swope 
telescope at Las Campanas Observatory, Chile. The camera uses a 2048 x 3150 SITe CCD 
with a pixel size of 15/im pixel"^ and a pixel scale of 0'.'435 pixel^^. A 1200 x 1200 section 
of the CCD corresponding to 8.7' on the side was read out, centered at a = 16:47:03.8, 5 = 
—45:50:41, J2000.0. The observations were made on 17 nights between UT 2006 June 15 
and July 25, typically in IRVB sequences. The total number of exposures were 252x5 s 
and 216x30 s in the / filter, 217x30 s in R, 200x600 s in V, and 167x1200 s in B. Two 
exposure times were taken in / to increase the sensitivity to both bright and faint cluster 
members. The median value of the seeing was 1.2" in /, 1.3" in i?, and 1.4" in both V and 
B. Westerlund 1 was observed at airmasses ranging from 1.04 to 2.69 in /, with the median 
being 1.16. 
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3. Data Reduction, Calibration and Astrometry 

The images were processed with standard IRAF0 routines. Specifically, the images 
were overscan corrected and trin imed, corrected for the non-linearity of the CCD with 
irred.irlincor in IRAF following iHamuy et al.l (120061 ) . and then fiat fielded. The value 
of 23,000 ADU, above which the non-linearity has not been measured, was adopted to be 
the saturation limit in the following reductions. 



3.1. Image Subtraction 



The phot ometry of the variable stars was extracted using the ISIS image subtrac- 
tion package (lAlard fc Luptonl Il998l : lAlardll2000l ). which works well in crowded fields (see 



Bonanos fc Stanek 20031) . In each filter, all the frames were transformed to a common co- 



ordinate grid and a reference image was created by stacking several frames with the best 
seeing. For each frame, the reference image was convolved with a kernel to match its PSF 
and then subtracted. On the subtracted images, the constant stars cancel out, and only 
the signal from variable stars remains. A median image of all the subtracted images was 
constructed, and the variable stars were identified as bright peaks on it by visual inspection. 
Finally, profile photometry was extracted from the subtracted images. 



DAOPHOT/ALLSTAR PSF photometry (jStetsoru 119871 ) was performed on the single 
template images sepa rately and the flux sca ling was corrected when transforming to magni- 
tudes, as described in lHartman et al.l (l2004j ) . A more detailed description of the calibration 
is presented in §3.2[ 



3.2. Photometric Calibration 



On two photometric nights, iLandoltl ( 1l992l ) standard fields were observed to calibrate 
the BVRI light curves to standard Johnson-Kron-Cousins photometric bands. Specifically, 
the SA107, SA 109, PG1323, PG1525 fields were observed on UT 2006 June 19 (N04), 
covering a range in airmass from 1.07 to 1.86, and the MarkA, PG0231, PG1633 fields on 
UT 2006 July 19 (N12), covering a range in airmass from 1.07 to 2.10. 

The transformation from the instrumental to the standard system was derived using the 
IRAF photcal package according to the equations: 



^IRAF is distributed by the National Optical Astronomy Observatory, which are operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under cooperative agreement with the NSF. 
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6 = 5 + Xb + «:b-X + eb-(B-V) 
Vi = V + Xvii + fi:vi2 ■ X + ^,13 ■ (B - V) 

V2 = V + Xv21 + fi:v22 ■ X + ^v23 " (V - I) 
r = + Xr + /tr ■ X + ■ (V - R) 

^ = / + Xi + «:i-X + ei-(V-I) 

where lowercase letters correspond to the instrumental magnitudes, uppercase letters to 
standard magnitudes, X is the airmass, x is the zeropoint, ^ the color and n the airmass 



coeffic ient. Additional bright red stars in the observed fields from the catalog of [Stetson 



(I2OOOI ) increased the color range for N04 to —0.22 < B — V < 1.90 mag using 95 stars and 
in N12 to —0.33 < B — V < 1.53 mag using 58 stars. Table [T] presents both photometric 
solutions. The solution of N04 has a larger number of stars, a greater color range and 
is in agreement wit h the values of the 5, V color and extinction coefficients measured by 



Hamuy et al.l (120061 ): it was therefore adopted. The transformation was repeated iteratively 
to determine the colors of the stars. The lack of observed standard stars with B — V > 
2 mag limits the accuracy of the photometry for the red dest stars. A com parison with 



the photometry of P. Stetson (priv. comm.) presented by IClark et al.l (120051 ) shows mean 
differences (Bonanos— Stetson) for the 100 brightest stars in common of —0.03 ± 0.04 mag 
in S, -0.01 ± 0.06 mag in V, 0.08 ± 0.08 mag in R and 0.24 ± 0.10 mag in I. The large 
offsets in R and especially in I are due to the lack of very red standards. The photometry 
was therefore scaled by these offsets in these two bands. 

The completeness of the photometry starts to drop rapidly at about 18.5 mag in /, 20 
mag in R, 21.5 mag in V and 22.5 mag in B. The CCD saturates for stars brighter than 
10.75 mag in / for the short exposures and 12.3 mag in the long ones, 12.7 mag in i?, 15.0 
mag in V and 16.4 mag in B. The light curves of variables brighter than 12.3 mag in / were 
measured from the short exposure images. 

The conversion of the variable star light curves to magnitudes involves two steps. First, 
it requires a conversion to instrumental magnitudes that takes into account the aperture 
correction measured on the template frame, to achieve correct scaling of fluxes and, conse- 
quently, correct amplitudes for the variables. This requires matching the variables with stars 
having DAOPHOT PSF photometry. In the few cases without matches, the light curves are 
presented in flux units. The conversion to the standard system involves an iterative color 
determination, which in turn requires detection of the variables in multiple fllters. This was 
not the case for a large fraction of the stars; therefore, the following cluster colors B — V = 3.8 
mag, V — I = 5.0 mag, V — R = 2.8 mag were adopted for the conversion, which are average 
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colors for the Westerlund 1 stars of interest. Hence, the photometry of the variable field 
stars is systematically offset by as much as 0.2 mag. Using these colors, the zeropoints were 
computed and added to the lightcurves. The I and R light curves were subsequently scaled 
to Stetson's photometry. 

3.3. Astrometry 

Equatorial coordinates were determined for the / star list. The transformation from 
rectangular to equato rial coordinates was derived using 940 transformation stars with V <2{] 
from the USNO-Bl.O (IMonet et al.ll2003l ) catalog. The median difference between the catalog 
and the computed coordinates for over 700 transformation stars was 0."2. 



4. Variable Star Catalog 

The BVRI light curves \ yere searched separately for v ariability using the multiharmonic 
analysis of variance method (jSchwarzenberg-Czernyl Il996l ) and then merged. A total of 129 
variables were found in the 8.7' x 8.7' region centered on Westerlund 1, many of these being 
members of the Westerlund 1 cluster. The variables were classified by visual inspection 
into the following five categories: (detached) eclipsing binaries or (D)EBs for binaries with 
periods P > 1 day, W Ursa Majoris (W UMa) contact binaries for P < 1 day and amplitudes 
> 0.2 mag, 5 Scuti for short period (P < 0.2 day) variables with amplitudes ~ 0.1 mag, 
periodic variables (P > 0.2 day), and "Other" that includes the long period or non-periodic 
variables. Table [2 ] presents the variable star catalog. It lists the na me of the star give n by 
Westerlundl ( 1l987l ). if available, or the revised identification given by lClark et al.l (120051 ) . the 
RA and Dec coordinates for J2000.0, th e P, V , R, I inte nsity- averaged magnitudes, the 
spectral classification when available from lClark et al.l (120051 ) and for the WR variables from 
Crowther et al.l (120061 ). t he light curve classif ication, and period when applicable. Variables 
with X-ray detections by ISkinner et al.l (120061 ) have an "X" following their WR c lassification. 



Note, that the catalog includes photometry for some stars that are missing from lClark et al. 



( 120051 ) due to the gap between their two CCDs and their smaller field of view. Table [3] 
presents the light curves of the variables. The name of each star is based on its J2000.0 
equatorial coordinates, in the format: hhmmss.ss+ddmmss.s. In five cases, the light curve 
is available only in differential flux units. 

The classification results yield 13 eclipsing binaries, 4 of which, and possibly 5, belong to 
Westerlund 1, 3 field W UMa binaries, 19 field S Scuti variables, 13 other variables exhibiting 
periodic variability, and 81 variables that are either long period variables or non-periodic 
variables. Because of the degeneracy in distinguishing between W UMa and S Scuti stars 
from photometry alone, a question mark is given with these classifications. SX Phe stars are 
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also short period pulsating stars in the classical instability strip, but belong to Population II 
and are unlikely to be found in the galactic disk. 

Figure [T] presents the J— band reference image (0."9 seeing) with the positions of the 
cluster eclipsing binaries and variable Wolf-Rayet stars. The eclipsing binaries belonging to 
Westerlund 1 were selected from the phased light curves as reddened eclipsing binaries with 
periods greater than a day. The location of the reddened DEB and WR77aa (T) demonstrate 
the extent of the cluster. 



4.1. Eclipsing Binaries in Westerlund 1 

Figure [2] presents phased light curves of the 5 newly discovered eclipsing binaries in 
Westerlund 1, in order of decrea sing brightness in I. The 9.2 day binary (W13, or star 13 in 



the catalog of I Westerlundl 119871 ) is more than a magnitude brighter than the other ec l ipsing 



systems in / (see Table |2]) and thus an excellent massive star candidate. I Clark et al.l (120051 ) 
assigned a spectroscopic identification of "OB binary/blend" for W13, which is in agreement 
with the contact eclipsing binary light curve shape. The second brightest system (W36), 
with the 3.18 day period, has a light curve resembling a near-contact configuration that is 
very accurate, thus lends itself to accurate parameter determination. There is no published 
spectroscopy for W36, but its red color and location, together with W13, in the core of the 
cluster give strong evidence for cluster membership. 



The third system, WR77o (B), is a suspected binary from X-ray observations. ISkinner et al 



(120061 ) detected strong, hard X-ray emission with Chandra observations as well as low- 
amplitude variability in its X-ray light curve and concluded it is a likely binary. These data 
confirm their suspicion, as well as show it to be an eclipsing binary with a 3.51 day period. 
The eclipses give it additional value, as they allow future measurement of the radii and 
masses of the component stars of this WN7 system. The uneven out-of-eclipse light curve 
implies non-uniform surface brightesses of the component stars, whereas the uneven eclipse 
depths imply different temperatures. 

The fourth system, a 4.43 day detached eclipsing binary in an eccentric orbit, is the 
first main sequence object to be identified in Westerlund 1, as the detached configuration 
requires unevolved stars. It is located far from the cluster core, ~ 4' south of the cluster, yet 
its colors are as reddened (V — / = 5.5 mag) as the cluster stars, indicating membership in 
the cluster. Similarly, WR77aa (T), at an even greater projected distance from the cluster 
center, is a member of Westerlund 1 (see Figure [1]). Measuring the masses of the above 
binary component stars will additionally determine the extent of mass segregation for the 
most massive stars in this super star cluster. Finally, the fifth eclipsing binary system 
found in Westerlund 1 has an extremely reddened color of R — I = 3.7 mag, is very faint 
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[R = 19.9 mag), and shows small-amplitude variability with a period of 2.26 days. Cluster 
membership of this star needs to be confirmed. It is located near the cluster core, however 
the i?— band light curve is too noisy to check whether the variation in amplitude is the same 
as in /—band, which is the case for eclipsing systems. It could alternatively be an extremely 
reddened pulsating or eclipsing field star. 

It would be premature to draw any conclusions about the binarity fraction from this 
data alone, since the monitoring baseline was only sensitive to short period eclipsing binaries. 
Spectroscopic monitoring will reveal additional non-eclipsing systems and, along with X- 



ray and radio observations, will yi eld a more comp 



ete e stimate of the binary fraction in 



(120061 ) to be > 60% for the Wolf-Rayet 



Westerlund 1, already estimated by ICrowther et al 
stars. 

4.2. Wolf-Rayet Variables in Westerlund 1 

Half of the 24 known Wolf-Rayet variables show significant optical variability, ranging 
from 0.05-0.4 mag. Most notably, WR77o (B), a WN7o type, is a 3.51 day period eclipsing 
binary with an eclipse depth of 0.2 mag, discussed in § I4.1[ The other 11 Wolf-Rayet variable 
stars are: WR77sc (A), WR77sb (O), WR77r (D), WR77q (R), WR77p (E), WR77n (F), 
WR77k (L), WR77j (G), WR77i (M), WR77f (S), WR77aa (T). Figures [3] and H present 
light curves of the 12 variable WR stars. ISkinner et al.l (120061 ) have detected X-ray emission 
in 9 of these 12 variable Wolf-Rayet stars, which correspond to the WR stars with the largest 
variability amplitudes and are marked with an "X" in Figures [3] and |H WR77f (S), with a 
WNIO — llh/BO — lla+ spectral type, is an exception. It varies by ~0.2 mag, yet is not 
an X-ray source. WR77sc (A) sh ows a sharp 0. 1 5 ma g increase in brightness that roughly 
phases with a period of 7.63 days. ISkinner et al.l (120061 ) conclude that both WR77o (B) and 
WR77sc (A) are binaries from their high X-ray luminosity. Eclipses are not visible in the 
light curve of WR77sc (A), however the periodic variability is likely related to winds. 

The low-amplitude variability (0.05-0.1 mag) observed in the rest of the WR stars 
qualitati vely resembles the li ght curve behavior of WR 123 (WN8), monitored by the MOST 
satellite (ILefevre et al.ll2005l ) and is likely due to inhomogeneities in the winds. The present 
data are not sufficient for determining whether a periodic signal is present, simil ar to the 9.8 



hour s ignal found in WR 123. The cause of this periodicity is not well understood. iDorfi et al. 



(120061 ) propose that strange-mode pulsations in hydrogen r ich envelopes of stars with high 
luminosity to mass ratios are a possible explanation, whereas iTownsend fc MacDonaldl (120061 ) 
argue for a model of g-mode pulsations driven by a deeper iron opacity bump. Future 
monitoring of Westerlund 1, yielding simultanesouly light curves of (at least) 24 Wolf-Rayet 
stars, will help in resolving the nature of the variability. 



- 9 - 



4.3. Other Variables 

The other periodic variables include 7 field detached eclipsing binary systems, 2 of which 
2 are eccentric, and 1 semi-detached or contact binary shown in Figure [5l The non-detached 
field system, with a 6.862 day period, exhibits an uneven light curve that suggests it is a 
magnetically active RS CVn binary. There are also 9 field small-amplitude variables with 
P > 0.2 days (see Figure [6]), 19 field small- amplitude variables with P < 0.2 (see Figure [7]), 
and 4 semi-regular variables. These include WR77sc (A) (P = 7.63 days) and W6, an OB 
supergiant displaying a 2.2 day periodicity with a 0.2 mag amplitude. Many of the S Scuti 
variables show significant scatter in their light curves, indicative of multiple mode pulsations 
that are typical in such stars. 

A large fraction of the remaining 81 variables classified as "Other" have red colors, and 
can therefore help in identifying additional supergiants and other massive stars that are 
cluster members. Figure [8] shows light curves of 6 selected cluster members: the LBV, the 
B[e] star, and four blue supergiants. 

5. Color Magnitude Diagram 

Figure M presents the R vs. R — I color-magnitude diagram (CMD), indicating the 
positions of the variables and in particular the eclipsing binaries in the cluster. Field stars 
are mainly located on the blue main sequence, whereas the Westerlund 1 stars have reddened 
colors betw' e en 2.2 < R — I < 3.2 mag. The sparse second main-sequence identified by 



Clark et al.l (120051 ) is also observed. It contains several variables, including the possibly 
magnetically active RS CVn eclipsing binary (P = 6.862 days), which is consistent with this 
group of stars being reddened by an intervening absorber. The brightest cluster eclipsing 
binary (W13) is remarkably bright, similar in brightness to the OB supergiants. Although 
a magnitude fainter, the second brightest EB (W36) is also located near OB supergiants on 
the CMD. The cluster DEB has a similar magnitude to WR77o (B); its position on the CMD 
is consistent with it being a cluster member. Spectral classification is required to confirm 
membership of the faintest EB. From their position on the CMD, the 4 semi-regular variables 
are likely cluster members; 2 of these (W6 and W72) have early-type spectral classifications 
and the brightest corresponds to W53. 

There is definite contamination from the classical instability strip, as several 6 Scuti 
variables occupy the same location as the cluster members. The 1.3 day variable with a 
0.1 mag J— band amplitude at P — J = 2.7, P = 15.9 mag could be a subgiant FK Comae 
star. These rapidly rotatin g G-K stars are thought to form from coalesced W UMa binaries 



see 



Bopp &: Stencellll98ll ). Finally, the bright periodic variables with periods of 3.7 days 
(P — / = 1.55, P = 16.3 mag) and 3.715 days (P — / = 0.9, P = 14.6 mag) could alternatively 
be low-inclination contact binaries with 7.4 day periods. 
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6. Summciry 

This paper presents the first variabiUty study of the Westerlund 1 cluster and is the first 
of a series of papers in search of the most massive stars (> SOM©) in echpsing binaries in 
young massive galactic clusters. A total of 129 variable stars are presented, including cluster 
supergiant and Wolf-Rayet star light curves, 4 (plus 1 candidate) cluster eclipsing binaries 
suitable for follow-up spectroscopy, 8 field eclipsing binaries, 19 pulsating S Scuti variables, 
3 W UMa eclipsing binaries, 13 other periodic variables and 81 long period or non-periodic 
variables. A significant number of the long period or non-periodic variables have red colors, 
suggesting cluster membership. 

The discovery of 4 or 5 eclipsing binaries belonging to Westerlund 1 is significant, as it 
provides the means of measuring accurate masses and radii for these stars. The brightest 
of these has V = 17.554, R = 14.707, / = 12.044 mag, therefore its radial velocity curve is 
measurable with 4-6 meter class telescopes. Future determination of the Wolf-Rayet eclipsing 
binary parameters will be especially useful for understanding WN7 spectral type stars and 
probing their evolutionary state. The cluster detached eclipsing binary is the first main 
sequence object to be identified in Westerlund 1 and its location outside the cluster core 
may provide insights into cluster formation and evolution. 

A systematic search for the most massive stars in eclipsing binaries in young massive 
clusters is efficient in discovering massive candidates, as shown to be the case in Westerlund 1. 
Radial velocity surveys are necessary for discovering non-eclipsing systems, however, they 
require large amounts of time on 6-8 meter class telescopes to identify binary candidates 
and obtain radial velocity curves. Photometric monitoring can be done using 1-2 meter 
telescopes and can easily identify short-period massive eclipsing binary candidates for follow- 
up spectroscopy. Eclipsing binary systems are thus extremely valuable probes of the most 
massive stars and can yield fundamental parameters of the most massive stars in both young 
massive clusters and nearby galaxies. 

I am very grateful to Craig Hcinkc for bringing Westerlund 1 to my attention and thus 
motivating this variability study. Many thanks to Peter Stetson for providing his unpublished 
photometry of Westerlund 1, to Lucas Maori for help with the photometric cahbration, 
to Kris Stanek for useful discussions and comments on the manuscript. I would like to 
thank Stella Kafka, Phil Massey and Tony Moffat for a careful reading of the manuscript. 
I acknowledge research and travel support from the Carnegie Institution of Washington 
through a Vera Rubin Fellowship. 



Fig. 1. — /—band finderchart of the observed 8.7' x 8.7' field of Westerlund 1 (0."9 seeing). 
Nortli is up and East to the left. Cluster eclipsing binaries are marked with red circles and 
variable Wolf-Rayet stars are marked with blue squares. Note, the large projected distances 
of the cluster DEB and WR77aa (T). 
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Fig. 2. — Phased light curves of the 5 echpsing binary members of Westerhmd 1, in order of 
decreasing brightness in /—band. The (RA, Dec) coordinates, periods and filters are labelled. 
The 3.51 day system corresponds to WR77o (B). Cluster membership of the faintest system 
needs to be confirmed. 
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Fig. 7. — Phased /—band light curves of a selection of 6 bright field S Scuti field variables. 
The (RA, Dec) coordinates and periods are labelled. The scatter indicates multiple mode 
pulsations that are typical in such stars. 
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Fig. 8. — Light curves of 6 cluster variables cl assified as "Other" . The (RA, Dec) coordinates 
are followed by the spectral classification and IWesterlundl (119871 ) name. The filter is labelled 
in the upper left hand corner. 
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Fig. 9. — i? vs. R — I CMD for Westerlund 1 showing the location of the variable stars with 
light curves in these bands. Cluster and field eclipsing binaries are marked with blue stars 
and blue filled triangles respectively, S Scuti with open green squares, W UMa with yellow 
filled trianges, other long period and semi-regular variables with cyan filled and open circles, 
respectively. The long period or non-periodic "Other" variables are marked with magenta 
squares. Many variables, in particular for field stars, were not measured in these bands and 
are therefore missing from the plot. 
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Table 1. Transformation coefficients for the Swope CCD. 



Instr. Mag. X i 



b7vo4(B-V) 2.981 ±0.009 0.205 ± 0.008 -0.058 ± 0.005 

V7vo4(B-V) 2.892 ±0.009 0.128 ±0.007 0.056 ± 0.005 

V7vo4(V-I) 2.888 ±0.010 0.130 ±0.008 0.053 ±0.005 

rAro4(V-R) 2.766 ±0.021 0.075 ± 0.017 -0.051 ± 0.014 

iA-()i(V-I) 3.216 ± 0.009 0.055 ± 0.007 -0.012 ± 0.004 

bivi2(B-V) 3.017 ±0.016 0.218 ± 0.010 -0.041 ± 0.006 

Vivi2(B-V) 2.890 ±0.016 0.161 ±0.010 0.087 ±0.007 

Vivi2(V-I) 2.886 ±0.016 0.164 ±0.010 0.071 ± 0.006 

rjvi2(V-R) 2.789 ±0.016 0.069 ± 0.010 0.007 ±0.010 

ijvi2(V-I) 3.189 ±0.014 0.092 ± 0.008 -0.040 ± 0.005 



Table 2. Westerlund 1 Variable Star Catalog 



Name 


ra 


Dec 




{V) 


(R) 


{!) 


Spectr. 


Class. 


Period 




(2000.0) 


(2000.0) 


(mag) 


(mag) 


(mag) 


(mag) 


Class.. 




(days) 




16:47:26.07 


-45:50:34.3 




18.794 








Other 






16:47:25.67 


-45:54:12.9 


18.881 


17.333 








Other 






16:47:20.83 


-45:46:34.1 




15.719 








Other 






16:47:20.75 


-45:50:02.6 




16.961 








Other 






16:47:20.37 


-45:48:54.7 




16.523 








Other 






16:47:19.50 


-45:49:57.3 




17.605 








Other 






16:47:18.72 


-45:47:58.3 


19.316 


17.758 








Per 


0.449 




16:47:18.51 


-45:50:30.7 


19.227 


17.450 


16.646 


15.474 




DEB 


2.50 




16:47:18.17 


-45:50:29.3 




19.387 








DEB 


6.834 




16:47:16.80 


-45:51:20.8 


19.456 


17.650 


16.781 


15.596 




DEB 


3.493 




16:47:14.19 


-45:48:00.7 




18.134 








Other 






16:47:13.62 


-45:48:57.0 


19.557 


17.289 


16.020 


14.563 




EB 


6.862 




16:47:13.39 


-45:49:10.5 






15.960 


12.990 




Other 






16:47:12.92 


-45:47:54.4 


21.296 


19.397 


18.471 


17.280 




WUMa? 


0.3657 




16:47:11.60 


-45:49:22.4 






16.211 


13.314 




Other 






16:47:11.06 


-45:49:30.9 






17.729 


15.018 




6 Scuti? 


0.144 




16:47:11.04 


-45:51:36.5 




20.268 








DEB 


3.24 




16:47:09.75 


-45:51:24.5 






19.702 


16.694 




S Scuti? 


0.134 


70 


16:47:09.39 


-45:50:49.6 








11.352 


B3Ia 


Other 




71 


16:47:08.45 


-45:50:49.3 




17.531 


14.158 




B3Ia 


Other 




72,WR77sc(A) 


16:47:08.35 
16:47:08.24 

16:47:07.84 


-45:50:45.3 
-45:50:56.4 

-45:51:47.9 






16.750 
17.123 

16.940 


13.728 
13.710 

13.939 


WN7b,X 


Per 
Other 
Other 


7.63 


WR77sb(0) 


16:47:07.68 


-45:52:35.6 






16.886 


14.047 


WN6o,X 


Other 




243 


16:47:07.52 
16:47:07.29 


-45:52:29.0 
-45:50:24.2 


19.873 


15.730 


16.811 


14.346 


LBV 


Other 
Other 




74 


16:47:07.08 


-45:50:12.9 






15.835 


13.312 




Other 





Table 2 — Continued 



Name 


RA 


Dec 


{B) 


{V) 


(i?) 


{I) 


Spectr. 


Class. 


Period 




(2000.0) 


(2000.0) 


(mag) 


(mag) 


(mag) 


(mag) 


Class.. 




(days) 




1d:47:07.01 


/I C AC\ AC\ 1 

—45:49:40.1 






1 O TAO 

16.702 


1 O AOO 

13.938 




(Jtner 






1d:47:0d.d8 


/I C /I T OO C 

—45:47:38.5 




1 O OCA 

18.250 


1 A A1 A 

14.919 






(Jtner 




16a 


i /IT A/^ o 

16:47:06.63 


A f P" A /I O 1 

—45:50:42.1 




16.248 






A OT 1 

A21a+ 


other 




15 


16:47:06.61 


/( C CA OA C 

—45:50:29.5 




1 A OA1 

19.301 






UB binary / blend 


Other 




Id 


1 f /IT A/^ A G 

1d:47:Uo.4d 


AC C A O/^ A 

—45:50:26.0 




1 T C C /I 

17.554 


1 A TAT 

14.707 


1 O A /I /I 

lz.U44 


(JB binary / blend ! 


hjB 


A OA 

9.20 




If* Ai"! C\f* O A 

16:47:06.34 


AC AO O/^ O 

— 45:48:zo.z 




OA AOC 

20.486 


1 A OOA 

19.230 


1 T TO O 

17.733 




TTl TA 

W U Ma ! 


A OT1 

0.371 


Zoo 


1 C /IT r\c OA 

16:47:06.29 


AC Af\ OO T 

—45:49:23.7 




1 T TOA 

17.739 






r 51a+ 


Other 




WR77r(D) 


16:47:06.28 


-45:51:26.4 






18.935 


15.600 


WN7o,X 


Other 






i ^1 T A/^ OT 

16:47:06.27 


AC r 1 AO A 

—45:51:03.9 






17.150 


14.065 




Per 


5.2 


14c,WK77q(Kj 


1 O /IT A/^ 1 A 

lo:47:(Jo.lU 


AC CC\ OO /I 

— 45:5U:z2.4 






1 T C O 

17.652 


1 c ooo 


WJN5o,A 


Other 




241,WK77p(l!>j 


1 /IT A/? AT 

16:47:06.07 


/( C CO AO O 

—45:52:08.3 




1 O OTA 

18.379 


1 C CTC 

15.675 


1 O A /I A 

12.94U 


WL9,A 


(Jtner 




25 


If /IT AC OO 

16:47:05.82 


AC C A O O O 

—45:50:33.2 




1 O OA C 

18.205 


1 C O O 1 

15.331 


1 O C TC 

12.575 




(Jtner 






1 /i T A C TA 

16:47:05.79 


AC CI O O O 

—45:51:33.3 






1 A A AT 

19.907 


16.182 






2.26 


18 


1 /I T A C TO 

16:47:05.73 


AC CC\ C (\ A 

—45:50:50.4 




17.610 


14.898 


12.228 




(Jtner 




d9 


1 /I T AC OA 

16:47:05.39 


AC C 1 OA C 

—45:51:30.5 




16.322 








other 




WK77o(_h5j 


16:47:05.37 


—45:51:04.7 






17.588 


14.367 


WiN7o,A 


Cjti 


3.51 


239,WR77n(F) 


16:47:05.21 
16:47:05.12 


-45:52:24.8 
-45:51:10.1 




18.130 


17.621 


14.354 


WC9d,X 


Other 
Other 




36 


16:47:05.08 


-45:50:55.1 


22.989 


18.975 


15.920 


13.199 




EB 


3.18 




16:47:04.99 


-45:50:13.3 






17.515 


14.782 




5 Scuti? 


0.1762 


19 


16:47:04.87 


-45:50:59.3 




18.784 


15.366 


12.452 


O9.5Ia-B0.5Ia 


Other 






16:47:04.77 


-45:49:47.1 


18.568 


16.896 








5 Scuti? 


0.115 


8a 


16:47:04.77 


-45:50:24.8 




15.913 






F5Ia+ 


Other 






16:47:04.69 


-45:52:06.8 






16.899 


14.405 




5 Scuti? 


0.1798 


28 


16:47:04.63 


-45:50:38.4 




17.146 


14.303 




O9.5Ia-B0.5Ia 


Other 




20 


16:47:04.68 
16:47:04.55 


-45:51:23.9 
-45:50:08.5 






15.660 




<M6I 


Other 
Other 





Table 2 — Continued 



Name 


RA 


Dec 


{B) {V) {R) (/) 


Spcctr. 


Class. Period 




(2000.0) 


(2000.0) 


(mag) (mag) (mag) (mag) 


Class.. 


(days) 





16:47:04 


52 


—45 


51 


19 


4 






16 


556 












Other 


238 


16:47:04 


41 


—45 


52 


27 


5 






17 


820 










09 5Ta-R0 5Ta 


Other 




1 6-47-04 


21 


—45 


51 

O -L 


07 


I 










1 5 


457 


12 


307 


WNQh- X 

WIN C/Ll.jJ\. 




9 


16:47:04 


17 


—45 


50 


31 


I 






18 


056 


14 


541 








Other 


33 


16:47:04 


13 


—45 


50 


48 


4 


20 


090 






12 


803 






B5Ia+ 


Other 


60 


16:47:04 


13 


-45 


51 


52 


1 










16 


041 


13 


384 


09 51a-B0 51a 


Other 


WR77ifG") 

vv ±v( i } 


16:47:04 


03 


—45 


51 


25 


1 










17 


893 


14 


561 


fWN7o')X 


Other 


66,WR77i(M) 


16:47:03 


96 


-45 


51 


37 


6 










16 


906 


13 


949 


WC9d 


Other 




16:47:03 


96 


-45 


52 


12 













17 


262 


14 


797 




A Scuti'^ 


46 


16:47:03 


93 


—45 


51 


19 


6 










15 


537 


12 


499 




Other 




16:47:03 


75 


—45 


51 


12 


7 










17 


042 


13 


879 




Other 


7 


16:47:03 


63 


—45 


50 


14 


2 






15 


972 










B5Ia+ 


Other 


43a 


16:47:03 


56 


-45 


50 


57 


6 










15 


378 


12 


647 


09 51a-B0 51a 


Other 


10 


16:47:03 


35 


-45 


50 


34 


6 


22 


720 


18 


516 


15 


257 






09 5Ia-B0 5Ia 


Other 


4za 


10.4/ .Uo 


/o 




ou 


oz 


i 














1 n 


(04 


TDCT„ 1 ? 

jjoia-|- 


utuer 




16:47:03 


24 


-45 


53 


11 


1 










18 


545 


15 


847 




5 Scuti? 




16:47:03 


11 


-45 


51 


31 


1 






19 


167 


15 


929 


13 


165 




Other 


237 


16:47:03 


10 


-45 


52 


18 


8 






19 


008 


13 


634 






iM6I 


Other 


6 


16:47:03 


06 


-45 


50 


23 


8 






18 


814 


15 


928 


13 


313 


O9.5Ia-B0.5Ia 


Per 


5,WR77f(S) 


16:47:02 


96 


-45 


50 


19 


7 


21 


552 


17 


792 


15 


035 


12 


508 


WN10-llh/B0-lIa+ 


Other 




16:47:02 


71 


-45 


50 


23 


9 


19 


651 
















Other 




16:47:02 


69 


-45 


50 


50 


3 










15 


340 


12 


669 




Other 


23a 


16:47:02 


57 


-45 


51 


08 


9 






18 


316 


14 


986 


12 


081 


O9.5Ia-B0.5Ia 


Other 


48 


16:47:02 


46 


-45 


51 


24 


9 










16 


818 


13 


966 




Other 




16:47:02 


32 


-45 


51 


15 


2 










17 


536 


14 


902 




Other 


61a 


16:47:02 


30 


-45 


51 


41 


8 






17 


506 


14 


633 






O9.5Ia-B0.5Ia 


Other 


11 


16:47:02 


25 


-45 


50 


47 


1 






17 


508 


14 


599 


11 


898 


O9.5Ia-B0.5Ia 


Other 



0.1337 



0.1581 



2.20 



Table 2 — Continued 



Name 


RA 


Dec 


{B) 


{V) 


{R) 


{I) 


Spcctr. 


Class. 


Period 




(2000.0) 


(2000.0) 


(mag) 


(mag) 


(mag) 


(mag) 


Class.. 




(days) 


12a 


16:47:02.23 


-45:50:59.0 




17.575 


13.611 




A5Ia+ 


Other 




49 


16:47:01.91 


-45:50:31.6 






16.392 


13.638 




Other 




52 


16:47:01.85 


-45:51:29.4 




17.876 


14.745 


12.065 




Per 


6.7 




16:47:01.69 


-45:52:57.8 






15.393 


12.966 




Other 




78 


16:47:01.56 


-45:49:57.9 




17.413 


14.621 






Other 






16:47:01.48 


-45:51:16.1 






18.240 


15.553 




5 Scuti? 


0.1571 


4 


16:47:01.42 


-45:50:37.4 


18.781 








F2Ia+ 


Other 






16:47:01.44 


-45:52:35.0 






15.334 


12.989 




Other 




57a 


16:47:01.36 


-45:51:45.6 




16.877 






B3Ia 


Other 




21 


16:47:01.13 


-45:51:13.4 






15.619 


12.861 




Other 






16:47:00.88 


-45:51:20.4 






17.336 


14.340 




Other 






16:47:00.52 


-45:48:29.7 






17.981 


15.806 




Per 


1.089 




16:47:00.54 


-45:52:23.3 


17.905 










Other 




53 


16:47:00.39 


-45:51:31.7 




18.950 


15.901 


13.198 




Per 


1.30 




16:46:59.93 


-45:50:25.3 






17.206 


14.818 




5 Scuti? 


0.1703 




16:46:59.77 


-45:51:40.0 






19.264 


16.536 




Other 




2a 


16:46:59.71 


-45:50:51.2 




16.978 


14.287 


11.771 


O9.51a-B0.51a 


Other 






16:46:58.97 


-45:50:17.8 


19.188 


17.650 








Other 






16:46:58.78 


-45:54:31.9 




20.313 


17.377 


14.816 




DEB 


4.43 


55 


16:46:58.42 


-45:51:31.4 




18.030 


15.283 


12.758 


O9.5Ia-B0.5Ia 


Other 






16:46:57.71 


-45:53:20.0 






14.741 


12.318 




Other 






16:46:56.89 


-45:52:04.4 






17.141 


14.917 




Other 






16:46:54.45 


-45:53:30.0 






17.408 


14.975 




5 Scuti? 


0.1271 




16:46:53.44 


-45:53:00.3 




18.061 


16.598 


14.425 




Other 






16:46:52.67 


-45:53:57.1 




16.804 


16.274 


14.724 




Per 


3.7 




16:46:52.24 


-45:51:15.5 


17.941 


16.546 








Other 






16:46:51.70 


-45:49:07.9 




16.289 








5 Scuti? 


0.4852 



Table 2 — Continued 



Name 


RA 


Dec 


{B) 


{V) 


(i?) 




Spcctr. 


Class. 


Period 




(2000.0) 


(2000.0) 


(mag) 


(mag) 


(mag) 


(mag) 


Class.. 




(days) 




16:46:50.83 


-45:53:02.7 






14.568 


13.700 




Per 


3.715 




16:46:50.05 


-45:51:05.4 




18.165 








Per 


1.211 




16:46:49.98 


-45:54:17.7 


21.469 


19.730 


18.874 


17.752 




EB 


0.5467 




16:46:49.82 


-45:52:08.6 




18.825 








Per 


2.36 




16:46:49.65 


-45:52:52.9 




18.887 


16.589 


14.406 




5 Scuti? 


0.1440 




16:46:48.38 


-45:53:56.6 






18.324 


15.682 




5 Scuti? 


0.1381 




16:46:47.28 


-45:50:50.5 


18.071 










Per 


0.626 




16:46:47.07 


-45:52:22.8 






16.187 


13.872 




S Scuti? 


0.1980 


WR77aa(T) 


16:46:46.30 
16:46:46.08 


-45:47:58.2 
-45:48:44.3 




17.372 


15.516 


13.063 


WC9d 


Other 
Other 






16:46:45.84 


-45:51:08.1 




18.808 


16.548 


14.415 




S Scuti? 


0.1708 




16:46:43.88 


-45:53:48.7 




19.646 








DEB 


2.229 




16:46:43.38 


—45:51:59.1 












DhiB 


5.78 




16:46:42.88 


-45:50:39.9 


18.383 










Other 






16:46:42.32 


-45:51:50.4 




18.115 








Other 






16:46:42.10 


-45:47:01.7 




21.391 


19.826 


18.233 




WUMa? 


0.398 




16:46:41.74 


-45:47:13.8 




16.647 


14.673 


13.773 




6 Scuti? 


0.1642 




16:46:41.65 


-45:46:41.1 


17.994 










5 Scuti? 


0.1841 




16:46:40.03 


-45:47:50.4 












S Scuti? 


0.1650 




16:46:40.08 


-45:51:07.8 




18.285 








Per 


0.738 




16:46:39.82 


-45:47:46.4 


16.940 










S Scuti? 


0.1642 



-28- 



Table 3. Light Curves of Variables in Westerlund 1 



Name 




Filter 


HJD-2450000 


mag/flux 


^rnag / ^ flux 


16:47:26.07- 


-45: 


:50:34.3 


V 


3901 


.50711 


18.798 


0.013 


16:47:26.07- 


-45: 


;50:34.3 


V 


3901 


.53865 


18.776 


0.011 


16:47:26.07- 


-45: 


:50:34.3 


V 


3901 


.59942 


18.808 


0.018 


16:47:26.07- 


-45: 


:50:34.3 


V 


3901 


.62506 


18.779 


0.018 


16:47:26.07- 


-45: 


:50:34.3 


V 


3901 


.65933 


18.778 


0.019 


16:47:26.07- 


-45: 


:50:34.3 


V 


3901 


.70142 


18.776 


0.019 


16:47:26.07- 


-45: 


:50:34.3 


V 


3901 


.74311 


18.788 


0.020 


16:47:26.07- 


-45: 


:50:34.3 


V 


3901 


.77163 


18.770 


0.019 


16:47:26.07- 


-45: 


;50:34.3 


V 


3901 


.79786 


18.774 


0.023 



Note. — Tabic; 3 is available in its entirety in the electronic version of the 
Astronomical Journal. A portion is shown here for guidance regarding its form 
and content. 



